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:onducting a 377 ° 2 proper motion survey in the ~V ; 
cool white dwarf contribution to the Galactic dark matter. Using the 250 ° 2 for which we possess three epochs, 
and applying selection criteria designed to isolate halo-type objects, we find no candidates in a 5500 pc 3 effective 
volume for old, fast My = 17 white dwarfs. We check the detection efficiency by cross-matching our catalogue 



Abstract. We are conducting a 377 ° 2 proper motion survey in the ^V and I bands in order to determine the 



with Luyten's NLTT catalogue. The halo white dwarf contribution cannot exceed 5% (95% C.L.) for objects with 



Mv = 17 and 1 < V — I < 1.5. The same con clusion applies to a 14 Gyr halo composed of white dwarfs with 



hydrogen atmosphere, as modeled by Chabrici 
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1. Introduction 

The rotation of the Milky Way is believed to be sustained 
by a massive dark halo, at least beyond one disk scale 
length. Old, cool white dwarfs (WDs) are one of the pro- 
posed constituents of that dark matter. A number of con- 
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straints have been placed on the total amount of halo 
white dwarfs (HWDs). One is based on the light emit- 
ted by WD progenitors ( |Graff et al., 1999 ). Extragalactic 
metallicity measurements also place strong constraints 
through the metal and helium enrichment due to the 
HWDs (Freese, 2000), although mechanisms have been 
proposed to evade those constraints. 

Over the past three years a great excitement has come 
from proper motion surveys used to find fast WDs in the 



Solar neighbourhood. The surveys conducted by Luytcn 



(1979a) and Knox et al. (1999) are jointly sensitive to 
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1500 pc 3 for M v = 17.5 HWDs flFlynn ct al., 200l[ ). They 
find no halo WD (Liebcrt et al., 1988), although this is 
disputed in the case of Luytcrj 's survey, whose efficiency 
is still a matter of strong debate, 23 years after publica- 
tion. On the other hand, several collaborations have re- 
ported proper motion detection of fast, cool WDs some 
of them being subsequently confirmed spectroscopically 



(Ibata et al., 2000, Oppcnhcimcr et al., 2001a). Nelson 



dards and Ogle ( Paczynski et al., 1999| ) secondary stan- 
dards in the Baade window: V E = V- (0.32 ± 0.03) (F- 1) 
and 1 E = I + (0.06 ± 0.03)(F - I). We estimate the 
precision of the colour transformations to be 0.1 mag for 
V — I =0-2, the range of expected colours for HWDs. 
Comparison between our Ze photometry and I as mea- 
sured by Denis (Fouque et al., 2000) reveals a small 



0.06, while a comparison between 



et al. ( !002 ) also reported the discovery of new HWD can- our photometry and the Photometric GSC ( Lasker et al. 
didatcs in the HST Groth-Westphal strip, although using 1988) reveals no significant offset, with a 0.05 mag uncer 



only two epochs. Finally, Monet ct al. (2000] ) discovered 
new high velocity WDs but the small sample of new ob- 
jects confirms the high efficiency of Luyten (1979b| ), thus 
decreasing the maximum density of halo WDs. 

On the microlensing side, Lasserre et al. (2000| ) put a 
30% upper limit (95% C.L.) on the contribution of HWDs 
of any age - or any compact object of similar mass - to 
a standard halo, while Alcock et al. (200"o| ) interpret their 
13 microlensing candidates towards the Large Magellanic 
Cloud as 0.15~0.9M Q objects contributing 20t"% to the 
standard halo. 

The interesting controversy developed around the 
Oppcnhcimcr et al. (2001h| ) results addresses the question 



of the nature of their findings: thick disk or halo objects. It 
is in any case agreed that these relatively bright Mi w 14, 
young WDs contribute very little to the Galactic dark 
matter, a few percent at most. It has been suggested how- 
ever that these could be the bright tail of a fainter, older 



population (Hansen, 2001). 

Our proper motion survey began in 1996 using the 
EROS 2 cameras. It was designed to study the contri- 
bution of older, fainter WDs to the halo, and also ad- 
dresses the question raised by Hansen (200l[ ). Like other 
proper motion searches, ours takes advantage of the large 
expected velocity dispersion of halo objects, and it is 
not intended to address the question of the low veloc- 
ity, or brighter, WDs that most recent surveys probe 
( pppcnhcimcr ct al., 20014 |Majewski fc Siegel, 2002] ). 
Here we report our results based on the analysis of 
250 ° 2 that were observed during at least three epochs. 



These yield a sensitivity about three times that of [Luyten 
(1979b), depending on HWD colours and luminosities. 

2. The EROS 2 proper motion survey 

2.1. Instrument and observations 



The EROS 2 wide-field imager ( [Bauer fc de Kat, 1998|) 
was designed to search for microlensing effects towards 
the Galactic bulge and disk and the Magellanic Clouds. 
Its two 1 02 CCD cameras are mounted at the Cassegrain 
focus of the 1-m Marly telescope at La Silla (Chile), with 
a pixel size of 0"6. These two 8k x 4fc mosaics are il- 
luminated through a dichroic beam splitter, which, to- 
gether with the CCD efficiency, defines the bandpasses. 
The colour transformations between the Eros system 
(Ve, Ze) and the standard Johnson-Cousins (V,I) sys- 
tem are determined by observing Landolt (1992| ) stan- 



tainty. 

Proper motion observations were performed one 
to two hours per dark night, within 90 minutes of 
the meridian to minimize atmospheric refraction. 
Our limiting magnitudes are about V — 21.5 and 
7 = 20.5. We observed 190 02 in the Southern Galactic 
Hemisphere (— 79 < 6 ga i < — 48 °), in the following 
strips along the a coordinate, Ad = 1.4 deg wide: 
22hl6min< a < 3h44min at <5 = -44 45 min, 
23h31min< a < lh34min at 5 = -40° 09 min and 
22h24min< a < 3h28min at 6 = -38 ° 45 min, and 
187 02 in the Northern Hemisphere (41 ° < 6 ga i < 59 °) 
with 10h57min< a < 13h23min at 5 = -12° and 
10 h 57 min < a < 12 h 53 min at S = -4 36 min. 

Only those fields (250 ° 2 ) with three epochs separated 
by one-year intervals (two-year total baseline), in one or 
two bands, are taken into account in this letter. 

2.2. Proper motion catalogue 

The reduction software for source detection, classifica- 
tion and catalogue matching was written in the frame- 
work of the EROS Peida++ package ( Ansari, 1996 ). Since 
photon noise dominates the astrometric errors for most 
of the search volume, we use a simple two-dimensional 
Gaussian PSF-fitting algorithm to determine stellar po- 
sitions. A rough star/galaxy classification is performed 
to limit galaxy contamination, with cuts chosen so that 
few stars are misclassified. The catalogues of the three 
epochs are geometrically aligned, with the deepest one 
taken as reference, using a linear transformation fitted to 
the 40 brightest stars, and matched within a search radius 
corresponding to a 6"yr _1 proper motion. The RMS dis- 
tance between matched stars provides an upper limit to 
the total astrometric error, which is 25 mas (Ict) for bright 
objects, degrading to 150 mas for V = 21 or I = 20. 

In order to minimize contamination by spurious can- 
didates we require that all objects be observed at least 
three times, separated by one-year intervals. The proper 
motions are derived from linear motion fits. The \ 2 distri- 
butions deviate from the expected ones only for confidence 
levels lower than 0.5%, indicating the onset of systematic 
effects. The 4% of objects so affected, mostly poorly mea- 
sured galaxies, stars close to CCD defects or visual bina- 
ries, are removed from the analysis. 

3. Search for halo white dwarfs 
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3.1. Selection of halo candidates 

The choice of the criteria used to distinguish halo objects 
from the much more numerous disk and thick disk stars, 
and from extragalactic objects, is of prime importance. We 
apply two cuts: We first require that the halo candidates 
detected in the Ve band exhibit a reduced proper motion, 
H V = V + 5 log T77 ^ T + 5 = My + 5 log TT ^ pT - 3.38 > 
22.5, where [i is the proper motion and v± the tangential 
velocity, and that those detected in Xe have Hj > 21.5. 
This cut select candidates that are both fast and faint (see 
Fig§. 

We then require a proper motion [i > 0''7yr . This 
cut removes objects whose proper motion reality can- 
not be confirmed by examining our images. It affects 
our sensitivity, depending on the mean WD magni- 
tude: fimin = O'^yr^ 1 corresponds to a transverse veloc- 
ity of v± = ^f^- x 200 km s -1 at distance d. A typical 
200 km s _1 transverse velocity is expected for halo objects, 
and d — 60 pc is the typical distance of WDs in our effec- 
tive volume for My = 17 stars. In fact, a detailed simu- 
lation shows that this cut reduces our sensitivity by 20% 
(resp. 45%) for My = 18 (resp. 17) HWDs. 

The proper motion cut was set to \l > O'^yr" 1 as a 
result of the following optimization: on the one hand this 
value is low enough to ensure a high detection efficiency 
of HWDs and, on the other hand, it is sufficiently high 
to eliminate candidates which are likely to be spurious 
because of a poor measurement of their proper motion. 
Indeed, when the proper motion cut is lowered to 0"6 yr _1 , 
10 candidates pass the cut, but for those whose proper mo- 
tions were measured in both bands (6 among the 10 can- 
didates) we find that for all of them the measured values 
of their visible and red proper motion are incompatible. 

No candidates survive the two cuts: Hy > 22.5 
or Hj > 21.5, and fi > O'^yr" 1 . Had we required 
/i > 0"8yr _1 , our upper limits on the HWD contribution 
to the halo given below would move up by less than 1% 
of the standard dark halo. 



3.2. Comparison with \Luyten\ 's catalogue 



We cross match our proper motion catalogue with the 
NLTT catalogue ( |Luyten, 1979b| ), by searching for bright 
stars in our catalogue with proper motion higher than 
0"lyr _1 , within 1' of Luyten positions. It turns out 
that 85% of Luyten 's stars are recovered with compatible 
proper motion and proper motion angle, which demon- 
strates our high detection efficiency for stars similar to 



Luyten 's stars. The 15% loss is partly due to the expected 



loss due to defects in the CCDs (4%) and to pointing 
dispersion between epochs (4%). Additionally, our proper 
motion measurements for Luyten 's stars are significant at 
the 3cr level in each band, and are compatible at the 2a 
level. This gives us confidence in our proper motion ac- 
curacy, as most Luyten stars have proper motion much 
lower than our 0''7yr _1 halo selection cut. However, this 
only applies to bright stars, as there are no Luyten stars 



H < 6" yr 



Halo selection zone 




17 18 19 
V magnitude 

Fig. 1. Proper motion vs. V magnitude for all objects de- 
tected in the Ve band, whose proper motion significance 
is higher than 99.7% C.L. (big dots), and for a simulated 
isotropic standard halo of My — 17 WDs fully composed 
of 0.6 M© white dwarfs (small dots, statistics xlO). The 
selection cuts exclude the hatched area. Here 45% of the 
HWDs are lost due to the /i > 0"7yr _1 cut. 



close to our photometric detection limit. Unfortunately we 
can draw no conclusion regarding Luyten 's completeness 
based on this study, as we lack a reliable estimate of the 
contamination of our disk-like proper motion catalogue. 
This study will be reported in more detail in a forthcom- 
ing article. 



4. Discussion 

4.1. Halo model predictions 

In order to check our sensitivity to HWDs, we perform a 
full simulation of our observations. We measure the de- 
tection efficiency of a star on each frame by adding sim- 
ulated stars to real images. Using the Besangon model 
of the Galaxy (Robin et al., 200C), to which we add one 
HWD per field, we then construct a simulate catalogue 
for each of our fields and each of our images, and process 
these catalogues in the same way as our observations. The 
Besangon model is a consistent hydrodynamical descrip- 
tion of Galactic star counts and proper motion, based on 
initial mass function and stellar evolution. 

We use various halo models with different kinematics. 
We use isotropic velocity dispersions of <7id = 100, 130 or 
156kms _1 , or anisotropic, spheroid-like dispersions, with 
o"id = 85, 105 or 125 km s -1 along the V (Galactic rota- 
tion) direction. We suppose halo rotations of —50, or 
50kms _1 . The expected number of detections depend on 
these parameters at the 5 to 10% level. For a given kine- 
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matics, flattened halos correspond to higher local densi- 
ties, and thus to proportionately more detections. 

We consider HWDs of 16.5 < M v < 18 and 
—0.5 < V — I < 1.5, in steps of 0.5 mag. We also 
gen erate HWDs wi th hydrogen atmosphere according to 
the Chabricr (1999 ) luminosity functions for halo ages of 
14 and 15 Gyr. The predicted number of detections, in 
either the Ve or the 2e band, are reported in Table 0, for 
the standard isotropic, isothermal halo, which has a local 
density of 8.10 -3 M Q pc -3 . Statistical errors are 5% while 
systematic errors are estimated to be less than 25%, the 
quadrature sum of errors in calibration (15%), efficiencies 
(10%) and parameters of the halo kinematics (15%). 



Table 1. Number of detected WDs, in cither Ve or 
2e band, predicted for a standard, isotropic halo with 
a local HWD mass density of 8.10~ 3 M pc~ 3 , and for 
0.6 M© HWDs of various colours and magnitudes, over 
250 02 . The effective volume is given between parenthesis 
for My = 17 in units of 10 3 pc 3 . 




10 & 



17 17.5 
magnitude M v 

Fig. 2. Constraints on the WD contribution to the stan- 
dard halo, for various WD colours, at the 95% C.L. The 
error bars correspond to the statistical error in the simu- 
lations. We exclude the hatched area. The empty box cor- 
responds to the two W Ds detected by |lbata et al. (2000 ), 
the star the result by Flynn et al. (2001) (no error bar 
given) . 
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Regarding Oppcnhcimcr et al. (2001a), we do not 
exclude a 1 to 2% contribution by bright M/=13-15 
HWDs. We stress that our selection cut, fi > 0"7yr _1 , 
is more severe than theirs, while preserving most of the 
halo sensitivity for the faint magnitudes of interest here. 
Rcid et al. (200l|), [Reyle et al. (200l[) and |Koopmans~fe 



We have not taken binarity effects into account. 
Double stars of similar magnitude would appear brighter 
than individual WDs, so that objects at the brighter, 
slower end of the distribution could be lost by the reduced 
proper motion cut. Furthermore, the number of systems 
to which we are sensitive would be smaller than the actual 
number of stars; our constraints would have to be relaxed 
accordingly. 



4.2. Constraints on the halo WD fraction and age 



The constraints on the contribution of WD to the halo 
are obtained from the above predictions and the lack of 
halo candidates, and shown in Fig. ^. We find that a 
14 Gyr halo cannot be made of more than 5% 0.6 or 
1.0 M Q HWDs, at the 95% C.L., while a 15 Gyr halo can- 
not be made of more than 15% of 0.6 M HWDs, or 45% 
of l.OM HWDs. 



4.3. Comparison with other surveys 

Our results, based on a larger volume than the study by 



Flynn et al. (2001), confirm their conclusion. They are 
marginally compatible with results of [bata et al. (2000| ), 
depending on the exact WD colours and absolute magni- 
tudes. 



Blandford (2002) also discuss this survey. We do constrain 
the proposition of ([Hansen, 2001 ) that the WDs detected 
by Oppcnhcimcr ct al. (2001a) be the bright tail of a 



fainter population. If this population has a halo-like kine- 
matics, then the older part of the luminosity function, 
older by several Gyrs, cannot make up more than « 7% 
of the halo. 

Finally, our constraints are fully compatible with 
the EROS microlensing results. The constraints from the 
present study are even stronger than those from microlens- 
ing, for HWDs of absolute magnitude Mv < 18 and also 
if the halo age is smaller than 15 Gyr (using the Chabriei 
(1999|) luminosity function for DAs). 



5. Conclusion 

In this letter we have presented the analysis of 250 ° 2 of 
the EROS proper motion survey, dedicated to the search 
of old, faint white dwarfs. We found no halo white dwarf 
candidates. We exclude a white dwarf contribution to the 
halo above 5% for M v = 17, 1 < V- I < 1.5, at the 
95% C.L. Assuming a |Chabrier (1999| ) cooling model, a 
14 Gyr halo cannot be made of more than 5% of WDs 
with hydrogen atmosphere. Our limits degrade to 20% for 
a 15 Gyr halo of IMq white dwarfs or for My = 18, 
< V - I < 1.5 white dwarfs. 
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